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Abstract. In this article, we study the (d — l)-volume and the covering numbers of 
the medial axis of a compact subset of R d . In general, this volume is infinite; however, 
the (d — l)-volume and covering numbers of a filtered medial axis (the ^t-medial axis) 
that is at distance greater than e from the compact set can be explicitely bounded. 
The behaviour of the bound we obtain with respect to fi, e and the covering numbers 
of K is optimal. 

From this result we deduce that the projection function on a compact subset K of R d 
depends continuously on the compact set K, in the L 1 sense. This implies in particular 
that Federer's curvature measure of a compact subset of R d with positive reach can 
be reliably estimated from a Hausdorff approximation of this subset, regardless of any 
regularity assumption on the approximating subset. 



1. Introduction 

We are interested in the following question: given a compact set K with positive 
reach, and a discrete approximation, is it possible to approximate Federer's curvature 
measures of (see [10] or §2.2 for a definition) knowing the discrete approximation only 
? A positive answer to this question has been given in [8] using convex analysis. In this 
article, we show that such a result can also be deduced from a careful study of the "size" 
— that is the covering numbers — of the medial axis. 

The notion of medial axis 1 has many applications in computer science. In image 
analysis and shape recognition, the skeleton of a shape is often used as an idealized 
version of the shape [18], that is known to have the same homotopy type as the original 
shape [15]. In the reconstruction of curves and surfaces from point cloud approximations, 
the distance to the medial axis provides a estimation of the size of the local features that 
can be used to give sampling conditions for provably correct reconstruction [1]. The 
flow associated with the distance function dx to a compact set K, that flows away from 
K toward local maxima of dx (that lie in the medial axis of K) can be used for shape 
segmentation [9]. The reader that is interested by the computation and stability of the 
medial axis with some of these applications in mind can refer to the survey [2]. 

The main technical ingredient needed for bounding the covering numbers of the subsets 
of the medial axis that we consider is a Lipschitz regularity result for the so-called normal 
distance to the medial axis. When K is a compact submanifold of class C 2 ' 1 , this function 
is globally Lipschitz on any r- level set of the distance function to K, when the radius r 
is small enough [12, 14, 5]. When K is the analytic boundary of a bounded domain 
of K 2 , the normal distance to the medial axis of dCl is 2/3-H61der on Q [3]. 



^also known as ambiguous locus in Riemannian geometry 
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However, without strong regularity assumption on the compact set K, it is hopeless 
to obtain a global Lipschitz regularity result for tk on a parallel set of K. Indeed, such a 
result would imply the finiteness of (d — l)-Hausdorff measure of the medial axis, which 
is known to be false — for instance, the medial axis of a generic compact set is dense. 

We show however, that the normal distance to the medial axis is Lipschitz on a suitable 
subset of a parallel set. This enables us to prove the following theorem on the covering 
numbers of the //-medial axis (see §3.1 for a definition): 

Theorem 4.1. For any compact set K C M. d , a parameter e smaller than the diameter 
of K , and rj small enough, 



From this theorem, we deduce a quantitative Hausdorff-stability results for projection 
function, which is the key to the stability of Federer's curvature measure (see Proposition 
2.2): 

Theorem 5.1. Let E be a bounded open set ofM. d . The application that maps a compact 
subset ofM. d to the projection function y>k £ L 1 (£') is locally h-Holder, for any exponent 
h smaller than 1 / (4<i — 2). 

Note that a similar result with a slightly better Holder exponent has been obtained 
in [8]. However, the proofs in this article give are very different and give a more geo- 
metric insight on the Hausdorff-stability of projection functions. Nonetheless, the main 
contribution of this article lies in Theorem 4.1. 



2.1. Distance, projection, boundary measures. Throughout this article, K will de- 
note a compact set in the Euclidean <i-space R rf , with no additional regularity assumption 
unless specified otherwise. The distance function to K, denoted by dx : X — > R + , is 
defined by dx{x) = m.m pe x \\p — x\\. A point p of K that realizes the minimum in the 
definition of dx{x) is called an orthogonal projection of x on K. The set of orthogonal 
projections of x on K is denoted by proj^-(x). 

The locus of the points which have more than one projection on K is called 

the medial axis of K. Denote this set by Med(K). For every point x of R d not lying in 
the medial axis of K, we let pk(%) be the unique orthogonal projection of x on K. This 
defines a map px ■ \ Med(K) — > K, which we will refer to as the projection function 
on the compact set K. 

Definition 2.1. Let K be a compact subset and £ be a measurable subset of R d . We 
will call boundary measure of K with respect to E the pushforward of the restriction 
of the Lebesgue measure to E on K by the projection function px, or more concisely 

VK,E = VK# ^ d \ E - 

We will be especially interested in the case where E is of the form K r , where K r 
denotes the r-tubular neig hborhood of K, i.e. K r = d^([0, r}). 




2. Boundary measures and medial axes 




Figure 1. Medial axis of a curve C in the plane, and Voronoi diagram 
of a point cloud P sampled on the curve. 

Example 2.1 (Steiner- Minkowski). If P is a convex solid polyhedron of M 3 , F its set 
of faces, E its set of edges and V its set of vertices, then the following formula holds: 

VP,pr = ^p + r^^f +r 2 Y J K(e)'H 1 \ e + r 3 Y J K ( v )Sv 

feF e£E v£V 

where K(e) is the angle between the normals of the faces adjacent to the edge e, and 
K(v) the solid angle formed by the normals of the faces adjacent to the vertex v. 

For a general convex polyhedra the measure nx,K r can similarly be written as a sum 
of weighted Hausdorff measures supported on the i-skeleton of K, whose local density is 
the local external dihedral angle. 

Example 2.2 (Weyl). Let M be a compact smooth hypersurface of M. d , and denote by 
<Jiip) is the ith elementary symmetric polynomial of the (d — 1) principal curvatures of 
M at a point p in M. Then, for any Borel subset B of and r small enough, the 
jix^ifr-measure of B can be written as 

d-l f . 

fJ-K,Kr(B) = V]const(i,d) / Ui{p)dM{p). 
i=0 Jbhm 

This formula can be generalized to submanifolds of any codimension [19]. 

2.2. Federer curvature measures and reach. Following Federer [10], we will call 
reach of a compact subset K of R. d the smallest distance between K and its medial axis, 
i.e. reach(X) = min xeMed ( A - ) d K {x). 

Generalizing Steiner-Minkowski and Weyl tubes formula, Federer proved that as long 
as r is smaller than the reach of K, the dependence in r of the boundary measure \ik,k t 
is a polynomial in r, of degree bounded by the ambient dimension d: 
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Theorem 2.1 (Federer). For any compact set K C R d with reach greater than R, there 
exists (d+ 1) uniquely defined (signed) measures &° K , . . . supported on K such that 
for any r ^ R, 

d 
i=0 

where uj^ is the volume of the k-dimensional unit sphere. 

These measures are uniquely defined and Federer calls them curvature measures of 
the compact set K. 

2.3. Stability of boundary and curvature measures. The question of the stability 
of boundary measures is a particular case of the more general question of geometric 
inference. Given a (discrete) approximation of a compact subset K of M. d , what amount 
of geometry and topology of K is it possible to recover ? In our case, the question is to 
bound the Wasserstein distance between the boundary measures of two compact subsets 
as a function of their Hausdorff distance. 

Recall that the Hausdorff distance between two compact subsets K and K' is de- 
fined by du(K,K') = \\dx — dx'W^- The Wasserstein distance (with exponent one) 
between two measures ^ and v with finite first moment on R d is defined by Wi (/x, v) = 
minx,y — Y \\] where the minimum is taken over all the couples of random variables 
X, Y whose law are /x and v respectively. 

Proposition 2.2. Let E be an open subset and K,K' be two compact subsets ofM. d . 
Then, 

w ( fJ-K,E Vk',e \ 1 I, || 

Proof. See [8, Proposition 3.1]. □ 

Hence, in order to obtain a Hausdorff stability result for boundary measures, one only 
needs to obtain a bound of the type \\px — Pk'IIl 1 ^) = °(^-h{K, K')). The possibility 
to estimate Federer's curvature measures from a discrete approximation can also be 
deduced from a L 1 stability result for projection functions (see [8, §4]). 

3. A FIRST NON-QUANTITATIVE STABILITY RESULT 

Intuitively, one expects that the projections pk{x) and pk'(x) of a point x on two 
Hausdorff-close compact subsets can differ dramatically only if x lie close to the medial 
axis of one of the compact sets. This makes it reasonable to expect a L 1 convergence 
property of the projections. However, since the medial axis of a compact subset of R d is 
generically dense (see [20] or [16, Proposition 1.2]), translating the above intuition into 
a proof isn't completely straightforward. 



SIZE OF THE MEDIAL AXIS AND STABILITY OF CURVATURE MEASURES 



5 



3.1. Semi-concavity of dx and /i-medial axis. The semi-concavity of the distance 
function to a compact set has been remarked and used in different contexts [11, 17, 4, 15]. 
More precisely, we will use the fact that for any compact subset K C M. d , the squared 
distance function to K is 1-concave. This is equivalent to the function vk '■ -> M,x H> 
|| x || 2 — d 2 K (x) being convex. Thanks to its semiconcavity one is able to define a notion of 
generalized gradient for the distance function dx, that is defined even at points where 
dx isn't differentiable. 

Given a compact set K C M d , the subdifferential of the distance function to K at 
a point x G W 1 is by definition the set of vectors v € M. d such that d 2 K {x + h) ^ 
d|^(x) + — A \\h\\ 2 for all h £R d . The subdifferential of d^ at a point x is denoted 
by d x dx, it is the convex hull of the set {{p — x)j \p — x\ ;p£ proj^(x)}. 

The gradient V x dx of the distance function dx at a point x 6 IR rf is defined as the 
vector of d x dx whose Euclidean norm is the smallest, or equivalently as the projection of 
the origin on d x dx (see [17] or [15]). Given a point x 6 M. d , denote by ^k{x) the center 
and ik{x) the radius of the smallest ball enclosing the set of orthogonal projections of 
m on K. Then, 

^7 A X-J K (x) 



d K {x) 

2 MN 1/2 
|V^|| = (l-^fS) = co ^) 



where is the (half) angle of the cone joining m to B(7^(m), r^(m)) 



3.2. /i-Medial axis of a compact set. The notion of [/,- medial axes and /x-critical 
point of the distance function to a compact subset K of M. d were introduced by Chazal, 
Cohen-Steiner and Lieutier in [6]. We recall the definitions and properties we will need 
later. 

A point x of M. d will be called a /i-critical point for the distance function to it" (with 
jj, ^ 0), or simply a /x-critical point of K if for every /i 6 R d , 

+ /i) < <&(x) + jU ||fc|| d K {x) + ||h|| 2 . 

The point x is /[/-critical iff the norm of the gradient || V x d.R-|| is at most K. The ^-medial 
axis Med^(K) of a compact set K C IR d is the set of //-critical points of the distance 
function. Is is easily seen that the medial axis is the union of all /U-medial axes, with 
/J, < 1: 

Med(K) = J Med M (^). 

Moreover, from the lower semicontinuity of the map x i— >■ ||V x dii-||, one obtains that for 
every // < 1, the /i- medial axis Med^iT) of K is a compact subset of R d . The main 
result of [6] that we will use is the following quantitative critical point stability theorem. 

Theorem 3.1 (Critical point stability theorem). Let K,K' be two compact sets with 
d}i(K,K') ^ e. For any point x in the fi-medial axis of K, there exists a point y in the 
fi' -medial axis of K' with // = fi + 2^ e/dx{x) and \\x — y\\ ^ 2-\J edxix). 
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3.3. A first non-quantitative stability result. The goal of this paragraph is to prove 
the following non-quantitative L 1 convergence result for projections: 



Proposition 3.2. // (K n ) Hausdorff converges to a compact K C W 1 , 
bounded open set E, lim^+oo \\pK n - p K \\ L ir E \ = 0. 



then for any 



In order to do so, for any L > 0, and two compact sets K and K' , we will denote 
Al{K, K') the set of points x of R d \ (K U K') whose projections on K and K' are at 
least at distance L, i.e. ||p/^(x) — p^'(x)|| ^ L. For technical reasons, we remove all 
points of the medial axes of K and K' from Al(K, K'). Since the Lebesgue measure 
both medial axes vanishes, this does not affect the measure of Al(K, K') . 

A consequence of the critical point stability theorem is that Al{K, K') lie close to the 
//-medial axis of K for a certain value of fi (this Lemma is similar to [7, Theorem 3.1]): 

Lemma 3.3. Let L > and K,K' be two compact sets and 5 ^ L/2 denote their Haus- 
dorff distance. Then for any positive radius R, one has 



A L (K, K') n K R C Med^K) 2 ^ 



with 



L — 5 
AR 




Proof. Let x be a point in Al(K,K') with dx(x) ^ R, and denote by p and p' its 
projections on K and K' respectively. By assumption, \\p — p'\\ is at least L. We let q 
be the projection of p' on the sphere S(x, d^(x)), and let Kq be the union of K and 
q. By hypothesis on the Hausdorff distance between K and K', the distance between p' 
and q is at most 5. Hence, du(K,Ko) is at most 25. 

By construction, the point x has two projections on Kq, and must belong to the Re- 
medial axis of Kq for some value of fi. Letting m be the midpoint of the segment \p, q], 
we are able to upper bound the value of fiQ-. 



Ho < 



iV^d^J 2 ^ cos 



m 



' i ii i 
/ ll x — P\ 



Since p, q belong to the sphere B(x, dx{x)), one has (p — q) _L (m — x) and \\x — p\\ 

1 1 2 1 1 1 1 1 2 

\\x — m\\ + j \\p — q\\ . This gives 



Ro ^ 




-1/2 



1 + 



L — 5 
2R 



-1/2 



To get the second inequality we used \\x — m\\ ^ R and \\p — q\\ ^ L — 5. 

In order to conclude, one only need to apply the critical point stability theorem 
(Theorem 3.1) to the compact sets K and Kq with dfj{K, Kq) ^ 25. Since x is in the 
//o-medial axis of Kq, there should exist a point y in Med^-fT) with \\x — y\\ ^ 2\fR5 
and ijl = Ro + ^yS/L. □ 
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Proof of Proposition 3.2. Fix L > 0, and suppose K and K' are given. One can decom- 
pose the set E between the set of points where the projections differ by at least L (i.e. 
Al(K, K') n E) and the remaining points. This gives the bound: 

\\VK> ~ P^IIli(b) < L7i d (E) + n d (A L (K, K') n E) diam(if U K') 

Now, take R = sup^ ||dx||, so that E is contained in the tubular neighborhood K R , 
and fix L = e/H d (E). Then, for 5 = d K (K,K') small enough (e.g. less than some <5o), 
the value of fi given in Lemma 3.3 is smaller than one. Denote by //o the value given by 
the lemma for 5q. Then 

\\PK'-PK\\ h i {E) ^e + H d {Med^ (K) 2 ^)di&m(KUK') (3.2) 

Being compact, Med^ (K) is the intersection of its tubular neighborhoods. Combining 
this with the outer-regularity of the Lebesgue measure gives: 

lim^Med^i-Q 2 ^) = 'H c! (Med Mo (K)) = 0. 
5— >0 

Putting this limit in equation (3.2) concludes the proof. □ 



4. Size and volume of the //-medial axis 

From the proof of Proposition 3.2, one can see that a way to get a quantitative stability 
of the projection functions is to control the volume of tubular neighborhoods of some 
part of the fi- medial axis. Recall that the e- covering number of a subset X C M. d is the 
minimum number N of points x±, . . . , xn such that X is contained in the union of balls 
U^ 1 B(xj,e). The following inequality is then straightforward: 

U d {X £ ) ^ H(B(0, e))M{X, e) (4.3) 

Our goal in this section is to obtain a bound on the covering numbers of the considered 
part of the /u- medial axis (see Theorem 4.1) that will allow to control the growth of the 
volume of its tubular neighborhoods. 

Because of its compactness, one could expect that the /i-medial axis of a well-behaved 
compact set will have finite % (i-1 -measure. This is not the case in general: if one 
considers a "comb", i.e. an infinite union of parallel segments of fixed length in R 2 , 
such as C = Uj e N*[0;l] x {2~*} Q M 2 (see Figure 2), the set of critical points of the 
distance fonction to C contains an imbricate comb. Hence % rf-1 (Med^(C)) is infinite for 
any /u > 0. 

However, for any positive e, the set of points of the //-medial axis of C that are e-away 
from C (that is Med tl (C)nM. d \C £ ) only contains a finite union of segments, and has finite 
% d ~ ^measure. The goal of this section is to prove (quantitatively) that this remains 
true for any compact set. Precisely, we have: 

Theorem 4.1. For any compact set K C M. d , e ^ diam(i^), and r\ small enough, 
M (Med^K) n (R d \ K £ ), 77) ^ Af(dK, e/2) O 



diam(K) 


d-i\ 


IvV 1 - 
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Figure 2. The "comb" and a part of its medial axis (dotted) 



In particular, one can bound the (d — 1) -volume of the fi-medial axis 



U d ~ x (Med^K) n (R d \ K £ )^j ^ M(dK,e/2) O 



diam(ET) 



-i d-r 



Remark (Sharpness of the bound). Let x,y be two points at distance D in M. d and 
K = {x,y}. Then, Med(-fT) is simply the medial hyperplane between x and y. A point 
m in Med(K) belongs to Med^{K) iff the cosine of the angle 9 = \^.{x — m,y — m) is 
at most /i. 

\\x — y\\ 2 diam(.ftT) 2 



cos 2 (0) 



d 2 A 



m) 



4d 2 f (m) 



Hence, cos(O) ^ /i iff dx(m-) ^ \ diam(K) / 'y / 1 — /i 2 . Let z denote the midpoint between 



x and y; then dx(m) 



\z — m\\ + diam(K) /4. Then, Med^(iT) is simply the inter- 



section of the ball centered at z and of radius diam(K) \J / (1 — /x 2 ) with the medial 
hyperplane. Hence, 



^ d - 1 (Med M (K)) = 



diam(i^T)/i 2 

7^7 



d-r 



This shows that the behaviour in diam(i^) and /i of the theorem is sharp as /x converges 
to one. 

4.1. Outline of the proof. In order to obtain the bound on the covering numbers of 
the the 2e-away /i-medial axis Med fl (K)n(M. d \K 2£ ) given in Theorem 4.1, we prove that 
this set can be written as the image of a part of the level set dK £ under the so-called 
normal projection on the medial axis £ :M. d \ K — > Med(K). 

The main difficulty is to obtain a Lipschitz regularity statement for the restriction 
of the map £ to & suitable subset of dK £ . There is no such statement for the whole 
surface dK £ in general. However, we are able to introduce a subset C dK £ whose 
image under £ cover the e-away /i-medial axis, and such that the restriction of £ to S e 
is Lipschitz. This is enough to conclude. 



4.2. Covering numbers of the /i-medial axis. We now proceed to the proof of The- 
orem 4.1. 
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Definition 4.1. For any point x G M. d , we define the normal distance of x to the medial 
axis as tr-(x) := inf{i ^ 0; x + tV x 6.k € Med(K)}. We will set tk(x) to zero at any 
point in K or in the medial axis Med (if). 

For any time t smaller than t(x), we denote by ^^(x) the point ^ l K (x) = x + tV x dx- 
Finally, for any x K, we let £k{x) be the first intersection of the half-ray starting at x 
with direction V^d^ with the medial axis. More precisely, we define £r{x) = ^ T ^f\x) G 
Med(K). 

Lemma 4.2. Ze£ m be a point of the medial axis Med(K) with d(x,K) > e, and x be a 
projection of m on dK £ . Then £(x) = m. 

Proof By definition of K £ , d(m, K) = d(m, K £ ) + e, so that the projection p of x on K 
must also be a projection of m on K. Hence, m, a; and p must be aligned. Since the open 
ball B(m,d(m,p)) does not intersect K, for any point y e]p, m[ the ball B(y,d(y,p)) 
intersects K only at p. In particular, by definition of the gradient, V x dx must be the 
unit vector directing ]p,m[, i.e. V ' x dx = (m — x)/d(m, x). Moreover, since [x,p[ is 
contained in the complement of the medial axis, t(x) must be equal to d(x,m). Finally 
one gets Vl/ T ( x )(x) = x + d(x, m)V x dx = m. □ 

This statement means in particular that 2e-away medial axis, that is Med(-fT) PI (M d \ 
K £ ), is contained in the image of the piece of hypersurface {x G dK £ ; tk{x) ^ e} by 
the map £. 

Recall that the radius of a set K C R d is the radius of the smallest ball enclosing 
K, while the diameter of K is the maximum distance between two points in K. The 
following inequality between the radius and the diameter is known as Jung's theorem 
[13]: radius (K) ^2(1 + 1/d) ^ diam(if). 

Lemma 4.3. For any point m in the (i-medial axis Med^(K), there exists two projections 
x,y 6 proj^(m) ofm on K such that the cosine of the angle \^-{x — m,y — m) is smaller 

than 



V 2 

Proof. We use the characterization of the gradient of the distance function given in 
equation (3.1). If B(7^(m), rjf(m)) denotes the smallest ball enclosing proj^(m), then 
fi 2 ^ 1 — r 2 K (m)/d 2 K (m). Using Jung's theorem and the definition of the diameter, there 
must exists two points x,y in proj^(m) whose distance r' is larger than \^2r x(m). The 
following bound on the cosine of the angle 9 = \^-{x — m,y — m) concludes the proof: 



□ 

Lemma 4.4. The maximum distance from a point in Med^(-ftT) to K is bounded by 

,1/2 



^diam(K)/(l- M 2 )' 



Proof. Let x,y be two orthogonal projections of m G Med M (iT) on K as given by the 
previous lemma. Then, using equation (4.4) , one obtains 

11-r- _ ,,ll 2 IA 

< (l + M 2 )/2. 



d 2 K {m) 
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Hence, & 2 K {m) ^ ^(1 — fj, 2 ) 1 \\x — y\\ 2 , which proves the result. □ 

Let us denote by the set of points x of the hypersurface dK e that satisfies the 
three conditions below: 

(i) the normal distance to the medial axis is bounded below: r(x) ^ e ; 

(ii) the image of x by I is in the jU-medial axis of K: £(x) £ Med^(-fT); 

(iii) there exists another projection y of m = £(x) on dK £ with 



cos I -Z(p — m, q — m) J ^ 



2 V 7 ^ V 2 

A reformulation of Lemmas 4.3 and 4.2 is the following corollary: 

Corollary 4.5. The image of by the map I covers the whole 2s-away fx-medial axis: 
£(S%) = Med^K) n (R d \ K 2e ) 

4.3. Lipschitz estimations for the map I. In this paragraph, we bound the Lipschitz 
constants of the restriction of the maps Vd^, r and (finally) I to the subset S £ C dK £ . 

First, let dK £,t be the set of points x in dK £ where the distance function is differen- 
tiable, and such that t(x) is bounded from below by t. In particular, notice that is 
contaiend in dK £ ' £ . The following Lemma proves that the functions VP* and V^d^ are 
Lipschitz on dK £,t : 

Lemma 4.6. (i) The restriction of VP* to dK £ ^ is (1 + t/e)- Lipschitz. 

(ii) The gradient of the distance function, x i-> V^d^, zs 3/e-Lipschitz on dK £ ' £ . 

Proof. Let x and x' be two points of dK £ with r(x), r(x') > t, p and p' their projections 
on if and y and y' their image by VP*. We let n = 1 + t/e be the scale factor between 
x — p and y — p, i.e.: 

(*) y' — U = u(x' — x) + (1 — — p) 
Using the fact that y projects to p, and the definition of u, we have: 

\\y - p\\ 2 ^ \\y -p || 2 = \\y - p\\ 2 + \\p - p\\ 2 + 2 {y - p\p - p) 

i.e.O ^ \\p — p'\\ 2 + 2u(x — p\p — p) 

le.{p - x\p - p) ^ -u" 1 \\p - p'\\ 2 

Summing this last inequality, the same inequality with primes and the equality {p' — 
p\p ~ p') = ~ \W ~ p\\ 2 gives 

(**) (x' — x\p — p) ^ (l — ti" 1 ) \\p' — p|| 2 
Using (*) and (**) we get the desired Lipschitz inequality 

|| y — y'|| 2 = u 2 ||x — x'|| 2 + (1 — u) 2 \\p — p\\ 2 + 2u(l — u)(x' — x\p — p) 
< u 2 \\x - x'f - (1 - u) 2 \\p' - pf < (1 + t/e) 2 \\x -x'f 

□ 

The second step is to prove that the restriction of r to the set is also Lipschitz. 
The technical core of the proof is contained in the following geometric lemma: 



SIZE OF THE MEDIAL AXIS AND STABILITY OF CURVATURE MEASURES 11 

Lemma 4.7. Let to denote the intersection time of the ray x$ + tv$ with the medial 
hyper-plane H Xo>yo between xq and another point yo> and t(x,v) the intersection time 
between the ray x + tv and H xyo . Then, assuming: 

Oi\\xQ-y Q \\ ^ (v Q \x Q -y Q ), (4.5) 

Ik-yolKA (4.6) 

||i7 — vq\\ ^ A \\x — xq\\ , (4-7) 

£^t(x ,y ) (4.8) 

one obtains the following bound: 

6 

t(x, v) ^ t(x , v ) H — 7 (1 + AD) \\x - xq\\ 
as soon as \\x — xq\\ is small enough (namely, smaller than ea 2 (l + 3AD) -1 ,). 

1 1 2 1 1 2 

Proof. We search the time t such that \\x + tv — x\\ = \\x + tv — j/o|| , i.e. 

4.2 II 1 1 2 || I ■ 2 , r>,; i \ . ,2 |i 1 1 2 

t \\v\\ =\\x — y \\ + 2t{x — yo\v) + t \\v\\ 

Hence, the intersection time is t(x,v) = \\x — 2/o|| 2 /2(2/o — x\v). The lower bound on 
t(xo,yo) translates as 

i II 1 1 2 -i 

1 x - yo\\ . 1 n n 

2 (x - yo\v Q ) 2a 

If V x t and V v t denote the gradients of this function in the direction of v and x, one 
has: 

./ \ 1 — y || 2 — ?/o) 
V v t(x,v; - 



V x t(x,v) 



2 (y - x\v) 2 
1 \\x - yofv + 2(yo - x\v)(x - 2/0) 



2 (2/0 - x|t>) 2 

Now, we bound the denominator of this expression: 

(x - y \v) = {x - y \v - v ) + (x - x \v ) + (x - yo\v ) 
^ a \\x - 2/0II - (1 + A \\x - 2/0II) \\x - x \\ 
^ a \\x — 2/0 1 1 — (2 + AD) \\x — xo\\ 
The scalar product (x — 2/0 1 v) will be larger than (say) ^ \\x — 2/0 1 1 provided that 

Qt 

(2 + AD) ||x - soil < 2 \\ x ~ 2/o II 

or, bounding from below ||x — 2/0II by ||xo — 2/0 1 1 ~~ l|xo — x\\ ^ 2ae — ||xo — x||, provided 
that: 

(3 + AD) ||x — xo|| ^ a 2 e 
This is the case in particular if ||x — xo|| ^ a 2 e(3 + AD)^ 1 . Under that assumption, 
we have the following bound on the norm of the gradient, from which the Lipschitz 
inequality follows: 

||Va.*(a;, < 6/a 2 and ||V„t(x,u)|| ^ 4D/a 2 

□ 
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Using this Lemma, we are able to show that the function t is locally Lipschitz on the 
subset Si C 8K £ : 

Proposition 4.8. The restriction of r to is locally L-Lipschitz, in the sense that if 
(x,y) £ are such that \\x — y\\ ^ 5o, then \\£(x) — £(y)\\ ^ L \\x — y\\ with 

L = ( ^" /£ ) ~'* = 0(«/i) 

In order to simplify the proof of this Proposition, we will make use of the following 
notation, where / is any function from X C M. d to R or R d : 

^Ps fix '■= sup{||/(z) - f(y)\\ I \\x - y\\ ; (x,y) £ X 2 and \\x - y\\ ^5}. 

Proof. We start the proof by evaluating the Lipschitz constant of the restriction of r to 
S^, using Lemma 4.7 (Step 1), and then deduce the Lipschitz estimate for the function 
t (Step 2). 

Step 1. Thanks to Lemma 4.3, for any x in S^, there exists another projection y of 
m = £(x) on dK e such that the cosine of the angle 9 = Z(x — m,y — m) is at most 
+ ^ 2 )/2- Let us denote by v = V^d^ the unit vector from x to m. The angle 
between yJc and v is ir/2 — 9. Then, 

, /l-u 2 \ 1/2 

cos(vr/2 - 9) = sin(0) = ^1 - cos 2 (0) ^ a := ( J 

As a consequence, with the a introduced above, one has a \\x — y\\ ^ a\(v\x — y)\. 
Moreover, ||x — y\\ is smaller than D = diam(i^ £ ) ^ di&m(K)+e. For any other point x' 
in S^, and v' = V x 'dx, one has \\v — v'\\ ^ A \\x — x'\\ with A = 3/e (thanks to Lemma 
4.6). 

These remarks allow us to apply Lemma 4.7. Using the notations of this lemma, one 
sees that t(x, v) is simply r(x) while t(x', v') is an upper bound for r(x'). This gives us: 

6 

r(x') T (x) + + XD) \\x - x'\\ 
^ t(x) + M\\x — x'\\ 
' \ + diam(K)/e\ 



where M = O 



as soon as x' is close enough to x. From the statement of Lemma 4.7, one sees that 
\\x — x'\\ ^ 5q with 5q = 0{e/M) is enough. Exchanging the role of x and x' , one proves 
that \t(x) — t(x')\ ^ M \\x — x'\\, provided that \\x — x'\\ ^ do. As a conclusion, 

Til n( 1 + diam(K)/e\ 
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Step 2. We can use the following decomposition of the difference £(x) — £{x'): 

£{x) - £(x') = (x' -x) + (t(x) - T(x'))V x d K + T(x')(V x d K - V x ,d K ) (4.10) 

in order to bound the (local) Lipschitz constant of the restriction of £ to from those 
computed earlier. One deduces from this equation that 



Lip 



5 



£U O + Lip, 



+ T 



Vd^| 



(4.11) 



Thanks to Lemma 4.4, one has |t(x)| = 0(diam(i"Q/(l — /i) 1 / 2 ); combining this with 
the estimate from Lemma 4.6 that Lip Vd^ls^ ^ 3/e, this gives 

Moo Lip 5o Wd K \J = 0(diam(^)/[e(l - p) 1 ' 2 }) (4.12) 
Putting the estimates (4.9) and (4.12) into (4.11) concludes the proof. □ 

In order to be able to deduce Theorem 4.1 from Proposition 4.8 we need the following 
bound on the covering numbers of a levelset dK r , where K is any compact set in M. d 
(see [8, Proposition 4.2]): 

M(dK r , e) ^ Af(dK, r)M(S d -\ e/2r) (4.13) 
Proof of Theorem 4-1- Applying Proposition 4.8, we get the existence of 



L = Lip, £\ SE = 0(diam(K)/(e v / l -//)) and 5 = 0{e/L) 
such that £ is locally L-Lipschitz. In particular, for any 77 smaller than Sq, 
M (Med M (K) n (R d \ K £ ), n )=N {£(S^, r,) 

^AT(dK £ / 2 , v /L). 

The bound on the covering number of the boundary of tubular neighborhoods (equation 
(4.13)) gives: 

(4.15) 



(4.14) 



N(dK £ / 2 , v /L) ^ N{dK,e/2)N (S d ~\ . 



yield 



V 



ILe 



d-l 



Equations (4.14) and (4.15), and the estimation N(S d 1 ,p) ~ ojd-ip" 

M (Ued^K) n {M. d \ K £ ),r^j = Af(dK,e/2) O 
Its suffices to replace L by its value from Proposition 4.8 to finish the proof. 

5. A QUANTITATIVE stability result for BOUNDARY MEASURES 



□ 



In this paragraph, we show how to use the bound on the covering numbers of the 
e-away /j-medial axis given in Theorem 4.1 in order to get a quantitative version of the 
L 1 convergence results for projections. Notice that the meaning of locally in the next 
statement could also be made quantitative using the same proof. 

Theorem 5.1. The map K i-> y>k £ L 1 (i?) is locally h-Holder for any exponent h < 
1 

2(2d-l) • 
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Proof. Remark first that if a point x is such that dx(x) ^ ^L—dn(K, K'), then by defini- 
tion of the Hausdorff distance, d^/(x) ^ \L. In particular, the orthogonal projections of 
x on K and K' are at distance at most L. Said otherwise, the set Al(K, K') is contained 
in the complementary of the ^ — 5 tubular neighborhood of K, with 8 := dn(A, K'). As 
in the previous proof, we will let R = ||dx|| £oo , so that E 1 is contained in the tubular 
neighborhood K R . 

We now choose L to be 5 h , where h > 0, and see for which values of /i we are able to 
get a converging bound. Using Lemma 3.3, we have: 



A L (K, K') HK R <Z (Med M (K) n (R d \ K^- 5 ^ 2 ^)^ 



For h < 1/2, the radius \{L — $) — 2\f~Rb will be greater than L/3 as soon as as soon as 
5 is small enough. Hence, 



A L (K, K') C\K R C (Med^K) n (K d \ i^ L/3 )) 



2VRS 



(5.16) 



The jU above, given by Lemma 3.3 can then be bounded as follows. Note that the 
constants in the "big O" will always positive in the remaining of the proof. 



M< 1 + 



L — 5 
AR 



-1/2 



+ 4VS/L 



= i + (-5 2h + <JV2-V2) 

The term will be asymptotically smaller than 1 provided that 2h < 1/2— h/2 i.e. h < 1/5, 
in which case \i = 1 — 0(5 2h ). By definition of the covering number, one has: 



H d (A L {K,K')nK R ) 



(Med^A^n (V\ A L / 3 )) 



2VRS 



N (Med^{K) n (V \ A L / 3 ) , 2yfR6) x 0{5 d ' 2 ) 



(5.17) 



The covering numbers of the intersection Med^(A) D (M. d \ A L / 3 ) can be bounded using 
Theorem 4.1: 



AA (Med M (A) n (V \ A L / 2 ) , 2\/fitf) 



= M(dK, L/2)0 



diam(Er)/v A R5 



(5.18) 



= AA(&A,L/2)0 (V^X ^ 1 )) 



Combining equations (5.17) and (5.18), and using the (crude) estimation M(dK, L/2) = 
0{l/L d ) = 0(5~ hd ), 

H d (A L (K, K') n K R ) < M(dK, L/2) o(<T fc (d-i)-*(d-i)+|d) 
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Hence, following the proof of Proposition 3.2, 

\\PK> ~ Pk\\ l i (e) ^ LU d {E) + n d (A L (K, K') n E) diam(K U K') 

The second term converges to zero as 5 = dn{K, K') does if h < 2 (2d-i) • ^ ms conciu des 
the proof. □ 
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